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I. INTRODUCTION 

Standard electromagnetic Compton scattering provides a unique tool for studying hadrons, which is one of the 
most fascinating frontiers of modern science. In this process the relevant Compton scattering amplitude probes the 
hadron structure by means of two quark electromagnetic currents. We argue that replacing one of the currents with 
the weak interaction current can promise a new insight. 

The paper is organized as follows. In SecHllwe briefly discuss the features of the handbag factorization scheme. We 
introduce a new set of phenomenological functions, known as generalized parton distributions (GPDs) [H, 0; S 3 01 , 
and discuss some of their basic properties in Sec. IIIII An application of the GPD formalism to the neutrino-induced 
deeply virtual Compton scattering in the kinematics relevant to future high-intensity neutrino experiments is given in 
Sec. HyI The cross section results are presented in Sec. |Vl Finally, in Sec. |Vl] we draw some conclusions and discuss 
future prospects. 

Some of the formal results in this paper have appeared in preliminary reports in Refs. [3] and whereas a compre- 
hensive analysis of the weak neutral and weak charged current DVCS reactions in collaboration with W. Melnitchouk 
and A. Radyushkin has been presented in Ref. Neutrino scattering off nucleons for neutral and charged currents 
was also discussed in Refs. flQl and fll!|, respectively. 

II. HANDBAG FACTORIZATION 

The underlying theory of strong interaction physics is the universally accepted nonAbelian gauge field theory 
known as quantum chromodynamics (QCD). The theory postulates that hadrons are composite objects made up 
of quarks, and that the color interaction between them is mediated by gluons as the gauge bosons. In principle, 
QCD embraces all phenomena of hadronic physics, in other words, the fundamental building blocks of hadrons (i.e. 
quarks and gluons) are known, and further the interactions between them are described by the well-established QCD 
Lagrangian. Nevertheless, knowing the first principles is not sufficient. The main difficulty is that QCD is formulated 
in terms of colored degrees-of- freedom, while the physical hadrons observed in the experiments are colorless. How to 
translate the information obtained from the experiments at the hadronic level into the language of quark and gluon 
fields has yet to be answered, and it represents a challenging task. 

One way to tackle this problem is to consider projecting these fields onto the hadronic states, and then use different 
phenomenological functions (e.g. form factors, parton distribution functions (PDFs), distribution amplitudes and 
GPDs) to describe the resulting hadronic matrix elements. The key idea of this approach is the handbag factorization. 
According to this property, in the specific kinematical regime featuring the presence of a large invariant, such as the 
virtuality of the probe or the momentum transfer, the amplitude of the process splits into the hard (short-distance) 
part with only one active parton propagating along the light-cone and the soft (long-distance) part represented 
by the lower blob, see Fig. [TJ The asymptotically free nature of QCD allows us to compute the short-distance 
interactions by means of perturbation theory, whereas the nonperturbative stages of interactions are expressed in 
terms of well-defined quark-gluon operators taken on the light-cone, and being sandwiched between the hadronic 
states. These matrix elements, which accumulate the information about the long-distance dynamics of the process, 
are process independent nonperturbative objects. Thus they are parametrized in terms of universal phenomenological 
functions, and they can be measured with the help of different probes (e.g. photons and weak interaction bosons). 

There are several situations, where the handbag contribution plays an essential role: 

• Both initial and final photons are highly virtual and have equal space-like virtualities. This situation corre- 
sponds to the forward virtual Compton scattering amplitude. Its imaginary part, through the optical theorem, 
determines the structure functions of deeply inelastic lepton scattering (DIS). 

• The condition on photon virtualities may be relaxed in a sense that the initial photon is still far off-shell but 
the final photon is real and the invariant momentum transfer to the hadron is small. The situation corresponds 
to the nonforward virtual Compton scattering amplitude, which is accessible through processes such as deeply 
virtual Compton scattering (DVCS) and deeply virtual meson production. 

• The configuration, in which both photons in the initial and final states are real but the invariant momentum 
transfer is large. The physical process corresponding to this situation is known as wide-angle real Compton 
scattering. 
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FIG. 1: Handbag diagram. In the lowest QCD approximation, both photons couple to the same quark line through point-like 
vertices. 



III. GENERALIZED PARTON DISTRIBUTIONS 

The description of specific aspects of the hadron structure is provided by several different phenomenological 
functions. The well-known examples are the so-called old phenomenological functions such as form factors, PDFs and 
distribution amplitudes since they have been around for a long time. On the other hand, the concept of GPDs is new. 
These new functions, as hybrids of the old ones, provide the most complete information about the hadron structure. 
In fact, the old phenomenological functions are just the limiting case of GPDs. Like PDFs, GPDs encapsulate 
longitudinal momentum distributions of hadron's constituents, however, at the same time, like form factors, they 
also provide the information about their transverse coordinate distributions, and hence give a comprehensive three- 
dimensional snapshot of the substructure of the hadron. Moreover, the universality of GPDs enables one to develop 
a unified description of wide variety of different hard (i.e. light-cone dominated) processes, both inclusive and exclusive. 

In the symmetric partonic picture, see Fig. O we treat both the initial and final hadron momenta symmetrically 
by introducing the average hadron momentum p = {pi + P2) /2. Here x is the usual light-cone momentum fraction. 
In addition, another scaHng variable is introduced, the so- called skewness ^. It is defined as the coefficient of 
the proportionality between the Hght-cone plus components of the momentum transfer, r = pi — p2, and average 
hadron momentum, ^ = r+/2p+, and it specifies the longitudinal momentum asymmetry. Clearly, in the skewed 
(nonforward) kinematics GPDs uncover much richer information about the hadron structure, which is not accessible 
in DIS, or in general in any inclusive process. By removing the parton with the light-cone momentum fraction a; -I- ^ 
from the hadron, and replacing it at some later point on the light-cone with the parton of fraction x — ^, GPDs 
measure the coherence between the two parton momentum states in the hadron as well as their spin correlations. 
Finally, by considering two different hadrons in the initial and final states one can study fiavor nondiagonal GPDs, 
for example, in the proton-to-neutron transition accessible through the exclusive charged pion electroproduction, the 
proton-to-A transition in the kaon electroproduction or the nucleon-to-A transition. 




FIG. 2: Symmetric partonic picture. 

Since —1 < x < 1, the momentum fractions a; ± ^ of the active partons can be either positive or negative, and for 
that reason GPDs have three distinct regions. When ^ < a; < 1 both partons represent quarks, while for — 1 < a; < — ^ 
they are both antiquarks. In these two regions GPDs are just a generalization of PDFs. In the central region, 
—i<x<£_, which is often referred to as the mesonic region, the partons represent the quark-antiquark pair. Here 
GPDs behave like meson distribution amplitudes. 

At the leading twist-2 level, the hadron structure information can be parametrized in terms of two unpolarized 
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and two polarized GPDs denoted by H, E and H, E, respectively. They are functions of three variables (a;,^,t), 
where t = is the invatiant momentum transfer, and further they are defined for each quark flavor /. GPDs have 
interesting properties linking them to PDFs and form factors. In the forward limit, pi — p2, r = Q, = 0, t = 0, they 
reduce to PDFs. In particular, Hf and Hf coincide with the quark density distribution /jv (x) and the quark helicity 
distribution Af^ (x). We write the so-called reduction formulas, 

^ ' ' ^ \ -In i-x) x<0 ^ ' 

and 

#n^,o,o)=( /y^/^\ (2) 

^ ^ ' ' ^ \ A/at i-x) X <0, ^ ' 

while Ef and Ef have no connections to PDFs. They are always accompanied with the momentum transfer, and 
hence simply not visible in DIS. Even though they have no analogue in the forward limit their limits do exist. In the 
local limit GPDs reduce to the standard nucleon elastic form factors (i.e. Dirac, Pauli, axial and pseudoscalar form 
factors) , 



J ^dx Hf{x,^,t)=Fif{t) , J ^dx Ef{x,^,t)^F, 



fit), (3) 



dx Hf {x, t) = gAf (t) , [ dx Ef (x, ^, t) = gpf {t) . (4) 



-1 J-i 
We call these relations the sum rules. 



GPDs in general are also relevant for the hadron spin structure. In terms of GPDs their second moments at i = 
give the quark total angular momentum, 

1 

^1 = 2^ dxx[Hf{x,^,t = 0) + Ef{x,^,t^O)]. (5) 

On the other hand, by decomposing Jq into the quark intrinsic spin AE (measured through polarized DIS) and quark 
orbital angular momentum Lq, Jq ~ AI]/2 + Lq, we can access Lq. Finally, since the total nucleon spin comes from 
quarks and gluons, 1/2 = + Jg, we can extract the total angular momentum Jg carried by gluons. 



IV. DEEPLY VIRTUAL NEUTRINO SCATTERING 



In recent years, significant effort was made to study GPDs through measurements of hard exclusive electroproduction 
processes. The simplest process in this respect is the standard electromagnetic DVCS process. Complementary to 
DVCS different combinations of quark flavors can be accessed by utilizing the weak interaction current in neutrino- 
induced DVCS (or, alternatively deeply virtual neutrino scattering - DVNS), namely 

iy{k)N{pi) e- {k')N' {j,2)i{q2) (6) 

for the charged current, and 

v{k)N{j,x) v{k')N{p2)i{q2) (7) 

for neutral current reactions. The presence of the axial part of the V — A interaction enables to probe a different set 
of GPDs, i.e. the charge conjugation odd combinations of GPDs as well as charge conjugation even combinations 
and hence independently measure both the valence and sea content of GPDs. Furthermore, one deals with a different 
flavor decomposition, in particular, less enhancement from u-quarks. An additional feature of the DVNS process 
is the study of GPDs that are not diagonal in quark flavors, such as those associated with the neutron-to-proton 
transition. 

In the most general case, there are three relevant diagrams for a DVCS-like process, see Fig. [31 The nucleon 
blob with the boson leg, which can be either a virtual photon or weak interaction boson, and photon leg represents 



5 



the virtual Compton scattering amplitude. This diagram is referred to as the DVCS diagram or the Compton 
contribution. The real photon can also be emitted by a lepton as depicted in the remaining two diagrams. They 
illustrate the Bethe-Heitler process, where the lower blob with one boson leg stands for the electroweak form factor, 
while the upper part in each diagram can be exactly calculable in QED. 




FIG. 3: DVCS (a) and Bethe-Heitler (b and c) diagrams contributing to the leptoproduction of a real photon. 



In the Bjorken regime, in which the virtuality of the weak boson —qf and the invariant mass of the weak 
boson-nucleon system, s = {pi + qi)'^ , are both sufficiently large while the Bjorken ratio Xb = —qf/[2{pi ■ qi)] is 
kept finite, the behavior of the virtual Compton scattering amplitude is dominated by the light-like distances. The 
dominant light-cone singularities are represented by two (s- and u-channel) handbag diagrams, see Fig. lU in which 
the hard quark propagator is convoluted with GPDs. In addition, by keeping the invariant momentum transfer 
to the nucleon as small as possible we arrive at the so-called DVCS kinematics, s > —qf ^ —t. Accordingly, the 
t-dependence enters only in the soft part of the amplitude. 




FIG. 4: The s- and M-channel handbag diagrams for the virtual Compton scattering amplitude. 



In our case in terms of symmetric momentum variables, p = {pi + P2) /2 and q = {qi + (72) /2, the virtual Compton 
scattering amplitude is given by the Fourier transform of the correlation function of the electromagnetic and weak 
(neutral or charged) interaction currents. 



w 



I / d*ze^''-^N'{p-r/2,s,)\T{J^^,{z/2)J^{~z/2)}\N{p + r/2,s,)} 



(8) 



The weak current couples to the quark current through two types of vertices, qqZ'^ and qqW^ , and hence the quark 
fields at coordinates ±z can either have the same or different fiavor quantum numbers. For that reason, we treat 
these two cases separately. The method employed to study the amplitude (HI in the DVCS kinematics is based on the 
nonlocal light-cone expansion of the time-ordered product of currents in terms of QCD bilocal operators in coordinate 
space 



^T{J^,,(z/2)J^J,{-z/2)} = --J-^^g/{4 [s'^'>-^O{-{z\0)-te>^''-^^Oi:^{z\0) 

^ ^ f 

- c{ [s^P'^'^Oi- (z |0) - ze^P''^Ol+ {z |0)] } , 



(9) 



where Qf denotes the electric charge of the quark with fiavor / (in units of |e|), Cy and are the corresponding weak 
vector and axial vector charges, s^'"^'^ = g^^P g^^ -\- g*^^^ gP^ — g^^^ gf^ is the symmetric tensor and e'*'"^'' the antisymmetric 
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tensor in Lorentz indices (/i, z^), and further 

Of |0) ^ [Vi/ (-2/2) i.V'/ (-^/2) ±{z^ -z)] , 
= [^/ {z/2) 7^75^/ {-z/2) ±{z^ -z)] 



(10) 
(11) 



are the vector and axial vector bilocal operators, respectively. Unlike the weak neutral current sector explicily 
presented here, the expansion for the weak charged and electromagnetic currents is a bit more complicated due to 
the fact that carries an electric charge. Thus we end up with different nucleons in the initial and final states and 
accordingly, the bilocal operators are accompanied with different electric charges as well as flavors. 

Next we isolate the twist-2 part of these operators, and sandwich it between the initial and final nucleon states. 
To construct the parametrization for nonforward matrix elements on the light-cone, we use a spectral representation 
in terms of eight GPDs, which have well-defined symmetry properties with respect to x. In particular, in case of the 
weak neutral current we have 

(iV(p2,S2)|0^'^(^|0)|7V(pi,si)),.=o = u{p2,S2) ^u{p^,s^) j' ^dx e'-^'^ Hf {x,^,t) 



(/ 



AM 



u{pi,si) / dxe''^P-'Ef{x,i,t) , (12) 



{N{p2,S2)\0(^{z\0)\N{pusi))^.^o = u{p2,S2) f(75u{pus,) ^dx e'^^'^ Hj {x,^,t) 



"(P2,S2)^75«(Pi,si) / dxe^^P-^Ej{x,tt) ■ 



(13) 



In contrast to the electromagnetic DVCS process, which measures only the plus GPDs (i.e. the sum of quark and 
antiquark distributions), the DVNS process gives also access to the minus distributions (i.e. the valence configuration) 
associated with the axial part of the V ~ A interaction. Moreover, it's important to note that since the matrix 
elements with the weak charged interaction current involve different fiavor combinations the corresponding GPDs, 
which parametrize these matrix elements, are nondiagonal in quark fiavor. 

The final expression for the weak neutral leading twist-2 Compton scattering amplitude reads 



4 (p • ?) 



(P • 92) 



{p^q^+p''qi^)-g^'' 



'^WN iC,t)u{p2,S2) li2U{pi,Si) +£^j^ {^,t)u{p2,S2) 



4M 



u(pi,si) 



-'^WN {i,t)u{p2,S2) ^275m(J31,Si) +£wN (^'0 ^ij^" (P2, S2) 75^ (Pl , Si) 



2M 



(p • 92) 



ie^-'^q^pPr, 



^WN (€' t) U {p2, S2) ^275M (Pl, Si) -£wN (?' 



(92 • r) 
2M 



■w(P2,S2)75W(pi,Sl) 



— - — - (^2 / — f 
-'^WN (^' *) " (P2, S2) ^2U (pi, Sl)-£y^,^ (^, t) U (P2, S2) ^ u (Pl, Si) 



(14) 



where 



H 



+(-) 

WN 



^WN^ (^' 



(x-^ + iO) f 
dx Hf {x, ^, t) 

dx 



{x,i,t) 
1 



± 



1 



(15) 
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- E t) (— ± -1^) , (16) 

- E ^/ t) T , (17) 

^^^^^ ^) - E Qf^viA) (,_7+zO) ^/" ^' 

= E ^?/4(.) £ t) (— T ^) , (18) 

are the so-called Compton form factors given by the integrals of GPDs. The amplitude lfT4|) has both real and 
imaginary parts. The real part is obtained by the principal value prescription, whereas the imaginary part constrains 
evaluating GPDs at the specific point x = 

In the laboratory frame as the target rest frame we choose a coordinate system, in which the weak boson four- 
momentum has no transverse components, see Fig.[5l The differential cross section is given by 



1 1 XBy"^ 1 1^,2 



dxBdQidtd^ 32 (27r)^ Q\ ^l + Ax^KP/Ql 



Tr , (19) 



where the invariant matrix element is the sum of both the Compton and Bethe-Heitler contributions, T = Tc ^T^bh- 
Furthermore, we require the following: 

• The energy of the incoming lepton beam is fixed at w = 20 GeV. 

• The invariant mass squared of the weak boson-nucleon system should be above the resonance region, s > 4 GeV^. 

• The virtuality of the incoming boson has to be large enough to ensure light-cone dominance of the scattering 
process, Q\ > 2.5 GeV^ 

• The momentum transfer squared should be as small as possible, e.g. 0.1 GeV^ < \t\ < 0.2 GeV^. 

We pick one kinematical point, in particular with the parameters uj — 20 GeV, Qf — 2.5 GeV^, xb — 0.35 and 
y — 0.19, and plot the invariant momentum transfer as a function of the angle Obj between the incoming virtual weak 
boson and outgoing real photon, see Fig. El The values vary from —0.15 GeV^ at zero degrees up to —1 GeV^ at 
0B-/ = 12°. Finally, we set the angle between the two planes to zero. In the so-called in-plane kinematics the polar 
angles of both incoming and scattered leptons are also fixed to ^ = 20.2° and 0' — 25.3°, respectively. 



8 




FIG. 5: The kinematics of the generalized DVCS process in the target rest frame. The lepton momenta k and k' define the 
lepton scattering plane and the momenta of the final photon and outgoing nucleon define the nucleon scattering plane. 
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FIG. 6: Invariant momentum transfer as a function of the scattering angle between the incoming virtual weak boson and 
outgoing real photon. 

We use a simple model for nucleon GPDs [H, [3, [IE, [l3| with the following assumptions: the sea quark 
contribution is neglibile and hence the plus and minus distributions coincide; the t-dependence of GPDs, which is 
driven by the corresponding nucleon form factors, factorizes from the dependence on the other two scaling variables, 
X and ^; the ^-dependence of GPDs appears only in the Ef distribution. 

In the weak neutral case we consider a neutrino- induced DVCS on a proton target through the exchange of Z°. 
There is no Bethe-Heitler contamination (the photon cannot be emitted from the neutrino leg) and for that reason, 
we measure the pure Compton contribution, 

= ^/2 |e| Gpu {k') 7, (1 - 75) u(fc)e; (92) T^'j, ■ (20) 

In the weak charged current sector we examine scattering off neutrinos, which convert into muons by emission of , 
on a neutron target. Here in addition to the Compton part. 



Ten = V2\e\GFu{k')^,{l-^^)u{k)el{q2)r{^ 



(21) 



9 



the Bethe-Heitler contribution is present, 

T^BHun = V2\e\GFe*,{q2)u{k') 

however, it's given only by diagram (b) of Fig. [3l Moreover, the isospin symmetry is used to express the nucleon 
matrix elements that are nondiagonal in quark flavor in terms of flavor diagonal matrix elements ■ 



ik' + q2f 



(22) 



V. RESULTS 



We present the cross section results for Ob^ < 12°, which corresponds to taking the invariant momentum transfer 
—t<—\ GeV^. The charged current Compton cross section is larger than for the neutral current interaction, 
see Fig. [71 We compare the magnitude of the Compton contribution relative to the corresponding Bethe-Heitler 
background by plotting both cross sections together on a logarithmic scale, see Fig. [H The Compton cross section 
is significantly larger than the Bethe-Heitler signal. Note also that for the Bethe-Heitler cross section one should 
expect a pole at the scattering angle Ob-^ = 25.3° when the scattered muon is collinear with outgoing real photon. 
On the other hand, in the electromagnetic current case the Bethe-Heitler contribution is larger than the Compton for 
0B-y > 3°, see Fig.[9l Finally, we plot the ratio of Compton versus Bethe-Heitler cross sections for both weak charged 
and electromagnetic DVCS. In the forward direction (i.e. for Os-y < 4°) the Bethe-Heitler cross section is strongly 
supressed by a factor more than 100 compared to the Compton contribution for the weak charged interaction current, 
whereas in the electromagnetic current case the ratio between the cross sections remains of order unity, see Fig. [Tol 
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FIG. 7: Compton cross sections for the weak neutral DVCS with the neutrino beam (solid line), electron beam (dashed line), 
and weak charged DVCS with the neutrino beam (dotted line). 
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FIG. 8: Compton (dashed line), Bethe-Heitler (dotted line), magnitude of interference (dashed-dotted line), and total (solid 
line) cross sections for weak charged DVCS. 
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FIG. 9: As in Fig.[8]but for electromagnetic DVCS. 
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FIG. 10: Ratio of Compton and Bethe-Heitler cross section for weak neutral DVCS with the electron beam (solid line), weak 
charged DVCS with the neutrino beam (dashed line), and electromagnetic DVCS (dotted line). 



VI. SUMMARY AND OUTLOOK 

GPDs provide the most complete and unified description of the internal quark-gluon structure of hadrons, which 
can be probed through wide variety of both inclusive and exclusive hard processes. 

DVNS is an important tool to complement the study of GPDs in the more familiar electro-induced DVCS. In 
addition to different combinations and flavor decomposition of GPDs, providing a direct measurement of their valence 
and sea contents, the process enables one to access the distributions that are nondiagonal in quark flavor. 

We have computed the twist-2 Compton scattering amplitudes for both weak neutral and weak charged interaction 
currents by means of the Hght-cone expansion of the current product in coordinate space. Using a simple model for 
the nucleon GPDs, which only includes the valence quark contribution, we gave prediction for cross sections in the 
kinematics relevant to future high-intensity neutrino experiments. 

UnHke the standard electromagnetic DVCS process, we flnd that at small scattering angles the Compton 
contribution is enhanced relative to the corresponding Bethe-Heitler contribution. Thus a contamination from the 
Bethe-Heitler background is less of a problem when extracting the weak DVCS signal. 

In the future we plan to use more reaHstic model for nucleon GPDs, which also includes sea quark effects, elaborate 
separately contributions from the plus and minus distributions, and include the twist-3 terms in order to apply the 
formalism at moderate energies. 
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